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Abstract 

We present a calculation of nuclear matter which goes beyond the usual quasi- 
particle approximation in that it includes part of the off-shell dependence of 

> 

the self-energy in the self-consistent solution of the single-particle spectrum. 

o 

^r) ' The spectral function is separated in contributions for energies above and be- 

O ' 

r*** ■ low the chemical potential. For holes we approximate the spectral function 

for energies below the chemical potential by a <5-function at the quasi-particle 
q ■ peak and retain the standard form for energies above the chemical potential. 

For particles a similar procedure is followed. The approximated spectral func- 
r^i ■ tion is consistently used at all levels of the calculation. Results for a model 

X 

calculation are presented, the main conclusion is that although several observ- 

ables are affected by the inclusion of the continuum contributions the physical 

consistency of the model does not improve with the improved self-consistency 

of the solution method. This in contrast to expectations based on the crucial 

role of self-consistency in the proofs of conservation laws. 
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I. INTRODUCTION 



Nuclear matter provides an ideal testing ground for the theory of strongly interacting 
many-body systems. Translational invariance leads to a particulary simple ground state 
structure where wave functions are given by plane waves. This allows to focus on the specific 
many-body scheme. Brueckner already identified the leading contributions by noting that 
the dynamics are dominated by contributions from ladder-type diagrams. Further progress 
was achieved by embedding of the equilibrium description of the ground state in a dynamic 
non-equilibrium theory M. This introduces the conservation laws of e.g. energy, particle 
number and momentum, even in a theory which only deals with static equilibrium where 
no explicit time development is present. Kadanoff and Baym showed that the requirement 
that the approximated theory preserves the conservation laws is not trivial 0. Closer 
examination reveals that the Brueckner approach is a subset of a 'conserving' approximation, 
the so-called T- matrix approximation p|. The first term of the Brueckner expansion ignores 
hole-hole propagation, which turn out to be crucial in preserving the conservation laws. 
Including these additional terms, as done e.g. in a non-relativistic approach in Ref. 
and in a relativistic model in Refs. ||, indeed shows a marked improvement in the physical 
consistency of the description of the many-body system and the preservation of conservation 
laws. However, deviations still persist and it is both necessary and interesting to investigate 
the causes. 

A crucial aspect in the proofs that the approximated theory preserves the conserva- 
tion laws is the self-consistency of the solution of the system. To make the solution of the 
T-matrix approximation feasible one has to resort to several additional practical approxima- 
tions. Note the essential difference between the approximation in the many-body scheme, 
as e.g. introduced by employing the T-matrix approximation, in contrast to the practical 
approximations needed to make a calculation feasible. One common simplification is to 
approximate the intermediate 2-body propagators by an angle-averaged form ||. Then one 
can decompose the T-matrix in partial waves and the spatial integration reduces to a one- 
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dimensional integration over the momentum variable and the full T-matrix is reconstructed 
by a sum over the separate partial wave contributions. As shown in Refs. H0] this is a 
reliable approximation. 

The integrations over the energy variable are usually simplified by assuming that the 
imaginary part of the self-energy is zero and taking for the real part the value of the self- 
energy at the pole of the propagator. The spectral function then reduces to a 5-function with 
its peak at the quasi-particle pole. This ensures that in this 'quasi-particle' approximation 
all energy integrations are easily carried out analytically. Together with the angle- averaging 
procedure the problem is then reduced to a self-consistent solution of the self-energy which 
now only depends on the magnitude of the momentum. 

In the self-consistent solution in the quasi-particle approximation of the T-matrix only 
the values of the self-energies at the quasi-particle peak are calculated. Clearly, a calcu- 
lation of the full energy-momentum dependence of the self-energy will check the validity 
of the quasi-particle approximation. If one finds a small imaginary part and therefore a 
weak dependence of the real part on the energy variable the quasi-particle approximation is 
probably justified. This was done at zero temperature in a non-relativistic approach by e.g. 
Ramos et al. || , in Ref . ||,[S| results for a relativistic Dirac-Brueckner model were presented. 
The Rostock group performed non-relativistic calculations at finite temperature, focussing 
on bound states and critical temperatures 0. All these calculations essentially find that 
most of the strength of the spectral function is indeed centered around the quasi-particle 
peak. However, a significant percentage (for nuclear matter in the range of 10-20%, depend- 
ing on the model) is moved to energies farther out. Since the equations are non-linear it 
is not a priori clear what the effects of incorporating this in the solution would be. In the 
present work we present an approximation of the spectral function that does include part 
of the off-shell dependence of the self-energy. We will investigate the effect of it on various 
observables in a practical calculation for a model interaction and compare the results with 
those of the quasi-particle approximation. The organisation of this paper is as follows. The 
next section is devoted to the theoretical framework, in the third section we present results 
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for a calculation using a model potential, while in the final section we will summarize our 
results and draw conclusions. 

II. THEORETICAL FRAMEWORK 

A complete account of the employed formalism and the philosophy behind it can be 
found in Ref. fH|. The central quantity in the many-body theory is the self-energy. It 
describes the effects of the medium and contains all the information on 1-body properties of 
the system. In uniform equilibrated nuclear matter it is a function of both the momentum 
and the energy. As any physically relevant function, it is an analytic function of the energy. 
Therefore the real and imaginary part of the self-energy are related by a dispersion relation. 
Since it is impossible to calculate the full theory, one has to resort to an approximation 
which sums a specific class of diagrams. In Refs. P|,PH| a method of obtaining systematic 
approximations for the self-energy is presented. One first establishes a Martin-Schwinger 
hierarchy of n-body Green's functions. These express lower-n Green's functions in higher-n 
ones. One then defines a formal self-energy in which all the information on correlations 
is present. The self-energy can be approximated in several ways, thereby truncating the 
Martin-Schwinger hierarchy and making the system solvable. Clearly, for the approximation 
of the self-energy constraints are set by the physical characteristics of the system. Several 
schemes that also respect conservation laws ('conserving approximations') have been dis- 
cussed in the literature flUJlOl . For nuclear matter with its strong short-range interaction the 
T-matrix approximation is appropriate. This is closely related to the Brueckner approach 
of summing ladder-type diagrams, but the T-matrix approximation goes beyond this by 
including also the propagation of hole-pairs. This is essential in fulfilling conservation laws 

Em- 

We now discuss to the representation of the T-matrix approximation in momentum 
space. We first introduce the propagators, using the Dyson equation we obtain the effective 
(dressed) advanced propagator in the well known form: 
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g {±) (uj,p) 



(2.1) 



w -£-E(±)(a;,p)' 

where the retarded self-energy is defined as EW(wj) = E(u; — ir),p),r] — » + . Inserting this 
in the definition of the spectral function one obtains the standard form 



a(v,p) = g (+) (uj,p) - g { \u } p) 



-Urn 



E<+>(w,p) 



(u; - £ - Jfe [E(+)(w,p)])a + (Jm (SW^p)])*' (2 ' 2) 
At temperature T the correlation propagators can be expressed in terms of the spectral 
function and the Fermi distribution f(oo) = (exp((a? — p)/T + l)^ 1 where \x is the chemical 
potential; 



g < {uj,p) = f{uj)a{uj,p) 
g > (uj,p) = (1 - f(u))a(u,p). 



(2.3) 



The set of propagators Eqs. |2~T| . |2~3| fully describes the 1-particle properties of the system. In 
equilibrium they are fully determined by the retarted self-energy £( + )(a;,p), reflecting again 
that the self-energy contains all the physical information of the system. In other words, in 
equilibrated nuclear matter only one independent function describes the system: specifying 
the spectral function determines the self-energy and vice-versa. 

In the T-matrix approximation the self-energy is given by (we use a notation where p is 
the four-momentum specified by the pair (u,p) and p denotes the usual 3-momentum 



S (+) (o;,p) 



^[(^(p-pO|t (+) (^ j p)I^(p-p , ))a^ < (^p / ) 

+(\(p-p')\T<(n,P)\±(p-p')) A g^(u J ',p') 



(2.4) 



where Q = to + u', P = p + p' . The T-matrices are given by: 

(p\T^\n,p)\p') = (p\v\p')- 

with the propagator 



|^ (P\v\f)g[V (n, P, p") (f l r« (n, ?) W) , (2.5) 
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doo 1 du; 2 g > (uj 1 , i+p)g > (u 2 , 1 -p) -g < {ui, ^+p)g < (uj 2 , f -p) 
(27r) 2 Q — oj\ — oj 2 + ze 



(2.6) 



and the correlation T matrix 




The set of coupled equations for the self-energy and the T-matrix are too complicated to 
solve in its full dependence on the spatial momentum and energy variable and one has to 
resort to the conventional approximations of angle- averaging the two-body propagator and 
the quasi-particle approximation discussed before. 

The equations above display the extensions of the T-matrix approximation over the 
Brueckner scheme in that it also propagates hole-pairs. This is also expressed in the equa- 
tions above. The Brueckner scheme is an expansion in terms of the density, in our formula- 
tion this is the term in first order of g K . Brueckner theory is retrieved by neglecting all higher 
order terms in g K to zero: the hole-hole propagating parts in g\ 2 (the terms with products of 
g < , they are of second order in g K ) and T < . Then g 12 only propagates particle-particle terms 
and T^ 4 ") only contains particle-particle ladder diagrams. Only the first term in the defini- 
tion of Eq. |2.4j contributes to the self-energy which is equal to the conventional Brueckner 
self-energy. A detailed account of the relation between the Green's function approach used 
in this work and Brueckner theory is found in Ref. 0. 



In the present study we go beyond the quasi-particle approximation and by including 
part of the off-shell energy dependence of the self-energy. In approximating the spectral 
function one has to be careful not to introduce unphysical aspects. One important property 
of the spectral function is the sum-rule: 



A. Sum-rule conserving approximations of the spectral function 




(2.8) 
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This just states that the spectral function is a probability distribution normalized to one. 
The interaction spreads the probability of finding a particle with a given momentum over 
the full energy range. The sum-rule ensures that the total strength is preserved. It is a 
straightforward consequence of the anti-commutation relations of the fields. However, it is 
also easy to prove that a sufficient condition for Eq. is that the underlying self-energy is 
analytic. If the imaginary part of the retarded self-energy is negative for all u in the upper 
half-plane, all poles of the retarded propagator will lie in the lower half-plane. From the 
definition of the spectral function it is clear that the sum of all the residues is then equal to 
the sum-rule. Using standard complex analysis this is equal to the residue at infinity, which 
is easily proven to be unity for physically reasonable self-energies. Hence, any physically 
meaningful approximation of the spectral function must conserve the sum-rule. If the sum- 
rule is violated (as e.g. in Ref. fLT|), the underlying self-energy is not analytic and thus not 
physical. Another good illustration of this point is the so-called Extended Quasi Particle 
approximation [|T2| , |i"3f . Essential in the proof that this approximation fulfills the sum- rule 
is the analyticity of the self-energy. 

Stability of the ground state at T = requires that the imaginary part of the self-energy 
at uj = \i vanishes 

Im[E (+) 0,p)] =0. (2.9) 

This results in a zero for the spectral function at u = /i. As was shown in Ref. || this 
property disappears for finite temperature, although remnants can still be seen at quite 
high energies. As we will see later, due to the analyticity of the self-energy this also has 
repercussions on the real part of the self-energy. At p = pf it results in the Fermi-surface 
consisting of ideal quasi-particles with infinite lifetime. This reflects the fact that the calcu- 
lated groundstate has to be stable: although all states in the Fermi sea have a finite lifetime, 
their decay is Pauli-blocked. This is directly related to the finding that the collision term 
in the equation of motion is zero in equilibrated matter at T = [TO] . Only the decay of 



particles at the Fermi surface is not Pauli blocked, but they have an infinite lifetime and the 



system is thus stable. 

The Extended Quasi Particle approximation (EQP) is an attractive candidate for a 
practical approximation of the spectral function that includes off-shell effects. Almost all of 
the energy integrations can be carried out analytically. Also the zero of the spectral function 
at uj — /i is preserved which is the case for more schematic approximations as e.g. using a 
single Lorentzian with the real and imaginary parts of the self-energy fixed to the values at 
the quasi-particle peak. However, as we will see later in a calculation, for a strong enough 
interaction the quasi-particle strength becomes larger than unity for low momenta. This is a 
consequence of the imaginary part of the self-energy being zero at u = /i and the self-energy 
being analytic. This feature makes the EQP difficult to implement in a nuclear medium. 

We seek an approximation that introduces off-shell effects as a perturbation on the 
quasi-particle solution. This is achieved by approximating the spectral function separately 
for particle and holes. We define particles as having momenta for which the quasi-particle 
peak lies above /x and holes have a quasi-particle peak below /i. In practice this amounts 
to splitting the spectral function into distinct parts for momenta higher or lower than the 
Fermi-momentum, respectively. We define a 'quasi-hole' by assinging the full hole-strength 
with energies smaller then \i into a single quasi-particle peak and retain the standard spectral 
function for energies above //. Since the spectral function is zero at uj — /i the approximated 
spectral function is still a continuous function of the energy. For particles we apply a similar 
procedure and we get 

a h (oj,p) = Z h 2-k5(uj - ujqp) + 9(uj - fx)a(u,p) 
a p (uj,p) = 9(/i — u)a(uj,p) + Z p 2tt5(u — u qp ) 

a(u,p) = 6(p f - p)a h (uj,p) + 0(p - p f )a p (uj,p). (2.10) 

The factors Z ph are determined by requiring that the approximated spectral function fulfills 
the sum-rule for all momenta. In the equations u qp denotes the position of the quasi-particle 
peak: the energy at which the real part of the denominator of the propagator vanishes. 
By construction the approximated spectral function conserves the same numerical values 
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for the momentum distribution n(p), defined by 

n {p) = / dua(u,p), (2-11) 

J — oo 

and the factors Z Pjh are simply expressed in terms of n(p): Z h = n(p), Z p = 1 — n(p). 

In the system and with the model interaction we consider in this work a few percent of 
the total strength is moved across /i, so the Z-factors are close to one and the quasi-particle 
peak still carries the bulk of the strength. This is exactly what we intended, namely to 
introduce the off-shell effects as a perturbation on the standard calculation in the quasi- 
particle approximation. 



B. Calculation of propagators and the T-matrices 

Inserting the approximated spectral function of Eq. |2.10| in the correlation propagators 
we obtain 

g < (uj,p) = 9(h-uj) [6(p f - p)Z h (p)2n5(uj - u qp ) + 9(p - p f )a(u,p)} 

g>(u, p) = -9{u - fi) [9(p - p f )Z p {p)2n5(uj - u qp ) + 9(p f - p)a{uj, p)} (2. 12) 



With these forms we can work out the two-body propagator in Eq. |2.5| . We find contributions 
of products of 5-functions, products of a 5-function with a tail of a spectral function and 
products of tails. For numerical simplicity we drop the latter contributions since they carry 
only a very small portion of the total strength (they are of order 0((1 — n(p)) 2 ) for holes). 
We have 

(+Vo p -\ - Z h(Pi) z h{P2)0(Pf ~ Pi)9(Pf - gg) ~ Z p {p x )Z p {p 2 )9{p x - p f )6{p 2 - p f ) 

il — oj 1 — u z 4- ie 

gp qp 



+ 
+ 



*» <fah a{uj 1 ,p 1 )Z h (p 2 )9(p 1 - Pf)8(pf - p 2 ) 
-oo 2tt Q — uj\ — Ug + ie 

M du 2 Z h (p 1 )a(u 2 ,p 2 )9(p f -Pi)d{p2 ~Pf) 
-oo 2tt tt — uj^ p — LU 2 + ie 

00 duh a(ui,pi)Z p (p 2 )9(pf - pi)9(p 2 -p f ) 
n 2ir — uji — uj^ p + ie 

00 du 2 Z p (p 1 )a(uj 2 ,p 2 )9(p 1 - Pf)9(p f - p 2 ) 



2ix Q- uji p - u 2 + ie 
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where we used the shorthand notation p 12 = f- ±p. On this propagator we apply a standard 
angle-averaging procedure ||. All factors including ^-functions in the momenta will reduce 
to the standard angle-averaged Pauli operators ||, Q pp , Qhh and Q p h = 1 — Q pp — Qhh- In a 
similar fashion also the Z-factors and the self-energies are replaced by angle-averaged values. 
The propagator then becomes 

(+Vn p -\ - z hh(P,p)Qhh(P,p) ~ Z pp (P,p)Q pp (P,p) 

9l2 \} l i^iP) — 



+ 



duj 2 Z h (pi)a(u;2,p2)Qph(P,p) 
-oo f il — uj x qp — uj 2 + ie 



OO 



du 2 Z p (p 1 )a(u 2 ,p 2 )Qph(P,p) 



n 7r Q - lo^ p - u 2 + ie 



(2.14) 



We will refer to the first part as the quasi-particle propagator, to the contributions from 
the integrations as the continuum contributions. The continuum contributions vanish for 
im(£) = 0, in that case all Z P) h factors are equal to unity and the quasi-particle propagator 
reduces to the conventional angle-average two-body propagator as e.g. used in Refs. [0,0. 
Angle-averaging the continuum contributions is achieved in the following way: for a given 
pair of (P,p) where Q p h(P,p) is non-zero the pair (pi,p 2 ) is formed for an intermediate 
value of the angle such that p± < Pf < p 2 . Then the integration over the energy variable 
is performed. Note that by construction the tails of the spectral functions appearing in the 
integrals are smooth functions. 

The imaginary part of the two-body propagator has now two contributions, one from 
the elastic cut in the quasi-particle propagator, the other from the imaginary parts of the 
continuum contributions. This is reflected by the unitarity relation obeyed by the T-matrix. 
From the T-matrix equation one derives the formal unitarity relation in partial-wave rep- 
resentation (in the case of coupled channels the summation over these in the intermediate 
states is implicit): 

im[(p\T^(n,p)\p')} = [ ^^(p|r( + )(fi,p)|/)/m^ 12 (fi,p,/)](/|r(-)(fi,p)b'). 

J TV 

(2.15) 
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Using Eq. |2.14| this relation can easily be worked out. The contribution from the quasi- 
particle part is trivial. Taking the imaginary part of the continuum contributions the de- 
nominator under the integral reduces to a 5 function which sets 102 = Q — uL, and we 
obtain 

Im[(p\T( + \n,P)\p')} = 

Po(Zhh(P,Po)Qhh(P,Po) ~ Z pp (P,p )Q pp (P,p )) , (+) , w ,^(_), n pv, a 
2^ + 2E'(P,p ) 9p {PIT } ^ P )\Po)(Po\T { \n,P)\p) 

r dn"r>" 2 r n 

- y -^-Qph{P,f)<n - ^ p)P2 ) [z h (p, Pl )e(/i - n + O - z p (p, Pl )6>(n - wj, - ^ 

x (p\T^\n,P)\p")(p"\T^\n,P)\p'), (2.16) 

where p = p denotes the momentum at the quasi-particle peak of gi%: Q + £- + ™ _|_ 
2E(P, po)gp = 0. £'(P, po)gp is the derivative with respect to the momentum of the self- 
energy: 

Z'(P,p ) qp = J^(uV,P,p)U. (2-17) 

At temperature T = and = 2/x, p has a value such that Qhh{P,Po) = Q PP (P,Po) — 

and the quasi-particle part does not generate an imaginary part. The same is true for the 

continuum part at O = 2/x. In the first ^-function, u) 1 < fi and the ^-function is zero. 

A similar argument holds for the second ^-function. Altogether, this results in Im[T(Q = 

2/i)] = 0. This leaves the possibility that the system develops a two-body bound state at 

this energy, which can only happen for Im(T) = 0. This was e.g. explored in Ref. ||. 

At T = the unitarity relation, Eq. |2.16| , provides a simple relation between T < and 
y(±) 

(p\T<(n,p)\p) = -e(2n-n)im[(p\T^\n,p)\p)], (2.18) 

and using Eq. [2.14| one finds 



( P \T<(n,p)\p) = 

Z m + Z H^P0) qp 
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r drl'rl' 2 

+ J J ^QAPy)a(tt-u 1 qp)P2 )Z h (P )Pl )6^-tt + u l qp ) 
x (p\T( + \n,P)\p")(p"\T^\n,P)\p>). (2.19) 

Again, we refer to the first part as the quasi-particle contribution and to the second as the 
continuum contribution. 



C. Evaluation of self-energies 

In calculating the self-energy we of course have to use the same approximated spectral 
function consistently in all expressions. In the first term of Eq. [2.4| , T^g < , this is straight- 
forward. For the second term the fact that Im\T < g^\uj 1 p)] = ^T < a(co,p) allows us to 
use the approximated spectral function. This provides us with the imaginary part of this 
contribution. Since both contributions are analytic functions we can then use a dispersion 
relation to calculate the real part of T < g <y ~\uj,p). This is facilitated by the fact that the 
imaginary part of this contribution is non-zero over a range of only a few hundred MeV 
around the Fermi energy. We can identify several contributions to the self-energy, originat- 
ing from various combinations of continuum and quasi-particle parts. Explicit expressions 
and more details of the calculation of the self-energy are given in the appendix. 

Now, we are almost in a position to calculate the self-energy. In the expressions given 
below we need T^ + > at all possible values of the variables (fl,P,p). This we achieve by a 
suitable interpolation. We separate the interpolation in one for the real part of T^ + \ and, 
using the unitarity relation the interpolation of the imaginary part is split into a quasi- 
particle part and a continuum contribution. In the interpolation of the quasi-particle part 
the factor with the Pauli-blockers in front is separated out. These trivial factors, which 
are easily calculated, depend strongly on the momenta and hinder a good interpolation 
numerically. It is more efficient to interpolate the term T^T^ and multiply with the 
factors later on, which greatly improves the stability of the results. 

After obtaining a self-consistent solution of the self-energy we can calculate the binding 
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energy and other observables. The energy of the system is given by 

1 r dtp f_ 

2 J (27r) 4V 2m Jy V ,PJ 

where / is the spin-isospin degeneracy factor, 1 = 2 for neutron matter and 1 = 4 for nuclear 
matter, respectively. The particle-density of the system is given by: 

^ 1 1 wr^' p) ^ 1 1 ^r n(p) - (2 - 21) 

In the quasi-particle approximation (and therefore also for a non-interacting system) this 
gives the standard relation between the Fermi momentum and the density: p = g^P/- A 
conserving approximation conserves the total particle number, and after switching on the 
interaction one must obtain the same result for the density. Deviations from this ideal situ- 
ation can have various causes. Excluding trivial numerical inaccuracies, deviations indicate 
insufficiencies of the approximation scheme. We stressed before that a crucial requirement 
for a conserving approximation is self-consistency. Although we improved upon this by in- 
cluding off-shell effects, our approach is still not fully self-consistent since we still ignore the 
width of the quasi-particle peak in our approximated spectral function. This might introduce 
deviations from the ideal particle conservation. On the other hand, in momentum space the 
conservation laws are only proven in the gradient expansion, which ignores gradients of the 
T-matrix in the Wigner transform ||10|| . So violations of the conservation laws also might 
point to inaccuracies introduced by the gradient expansion. Note however that we still can 
define an unique value of the Fermi-momentum by setting it to the momentum for which 
the quasi-particle energy is equal to the chemical potential. 



III. RESULTS 
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A. Details of the numerical calculation 



All the results presented in this study were calculated with a non-relativistic version 



of the Bonn-C potential [R], some numerical details are discussed in the appendix. This 
provides an adequate fit to the phase-shifts and leaves the possibility of exploring relativistic 
effects in the future. To keep the numerical effort within limits we have to restrict ourselves 
to the most important partial-wave contributions. In the isospin T = 1 channel we keep 
the 1 S'o partial wave, in the isospin T = channel the coupled-channel 3 S\ — 3 D\ partial 
wave. The tensor force of the latter channel has two effects. Firstly, it is very effective in 
moving strength of the spectral function across u = fi, thereby reducing n(p) and increasing 
the importance of the tail of the spectral function. On the other hand, however, it is also 
responsible for the presence of a bound state, discussed extensively in Refs. P,P, |l7^ , |T6 



Although a treatment of this 'deuteron-like' bound state is very well possible within our 
approximation scheme, we leave this for a future study. Since we are interested in the 
effects of the continuum contributions, we want to keep as much of the tensor force as 
possible. Reducing its strength by a factor 1/2 removes the bound state and we can perform 
a calculation in a range of densities around the nuclear matter saturation density. 

We first calculate a self-consistent spectrum within the quasi-particle approximation. 
Then, still in the quasi-particle approximation, the self-energies are calculated on a mesh in 
the energy-momentum plane, which, in turn, is used as a starting point for the full calculation 
including the continuum contributions. For nuclear matter it turns out that for an adequate 
fulfillment of the sum-rule Eq. |2.8| we need to include energies up to about 2 GeV. Using 
the change in the resulting quasi-particle spectrum as an iteration parameter we find rapid 
convergence: for depletions in the order of 6 % two iterations are usually sufficient. 
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B. Quasi-particle spectrum and self-energies 



The resulting quasi-particle spectra for a calculation at normal nuclear matter density 
(pf = 0.27 GeV/c) are displayed in Fig. |I[ Deep in the Fermi sea the real part of the self- 
energy in the continuum calculation is some 3 MeV less negative than the one calculated 
in the quasi-particle approximation. For larger momenta the difference gets smaller: at 
the Fermi momentum the difference is around 2 MeV, at twice the Fermi momentum it is 
negligible. For the imaginary parts we again find the largest difference for small momenta. 
The imaginary part of the self-energy in the continuum calculation is 1.5 MeV smaller (less 
negative) than in the quasi-particle approximation. Approaching the Fermi momentum the 
difference obviously vanishes since the quasi-particle energy at the Fermi surface equals the 
chemical potential resulting in a vanishing imaginary part of the self-energy. For higher 
momenta we do not find significant differences within the calculated range. 

In figure [| we show the full dependence on both the energy variable u and the momentum 
p of the real and imaginary parts of the self-consistent self-energy calculated at normal 
nuclear matter density. For a better presentation the projection angles are taken differently 
and the imaginary part is multiplied by a factor -1. At higher energies the imaginary part 
shows a similar behaviour for all momenta. The zero of the imaginary part at uj = fi is 
clearly visible. At lower momenta we observe a peak, this is due to the large phase shifts for 
low laboratory energies in the L = partial waves. The imaginary part can be viewed as a 
cross section averaged over the momenta. In the equation of motion these are multiplied by 
phase-space factors which ensure Pauli-blo eking. For low pi a b one predominantly averages 
over low total momenta, these are equivalent to low laboratory energies and one picks up 
the large cross sections. At higher pi a b the average is over higher momenta where the cross- 
sections are smaller. Together with the requirement that 7m[E(u; = //)] = one obtains 
the structure as observed in the figure. Since the self-energy is an analytic function and 
therefore observes a dispersion relation the structure seen in the imaginary part is reflected 
in the real part. This is clearly the case, with the largest structure being present at the 
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lower momenta. For momenta around the Fermi momentum the quasi-particle point lies 
at the left side of the ridge (from the point of view in the figure), for lower momenta it 
lies on the right side of the ridge. In the latter case, the derivative of the real part of the 
self-energy with respect to the energy is positive, leading to a quasi-particle strength larger 
than one. This is reflected by the spectral function having two maxima, centered around the 
quasi-particle pole. If the derivative gets larger then unity the spectral function even has 
two quasi-particle poles. In our system this can happen at lower densities, there the average 
of the momentum is taken over a smaller volume, leading to a larger imaginary part. Also 
the energy scale is compressed a bit (i.e. the maximum in the imaginary part lies closer to 
uj = fi). This all leads to a larger positive derivative at low momenta. As we already noted 
before this makes the application of the Extended Quasi-Particle approximation in a scheme 
similar to the one we employ problematic. A similar behaviour was found in Refs. p|,|T3 



and with the arguments given above we think this is an unavoidable feature of any realistic 
interaction having large cross-sections at low energies. 

In Fig. [3] we compare the energy dependence of the self-energy at a fixed momentum. 
We now use the variable uj' = u — u qp , defined as the distance to the quasi-particle energy. 
This allows for a meaningful comparison between the results, even when the underlying 
self-consistent quasi-particle spectra are different. The imaginary parts of the self-energies 
calculated in both approaches show only minor differences on the scale of the figure. Only 
the maximum below uj' = is slightly larger for the quasi-particle calculation. The difference 
between the real parts is nearly constant (as one would expect when the imaginary parts do 
not differ on the same scale). This is caused by a smaller imaginary part of the self-energy 
in the continuum calculation for higher energies. 

C. Binding energy and other observables 

Having found a self-consistent self-energy we can calculate the whole set of one-particle 
observables. In Tables [I] and |TJ we present values for several observables calculated in the 
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quasi-particle and in the continuum method of this paper, respectively. In the second column 
the binding energy is shown. For the continuum method these are in the order of 1 MeV 
higher than the ones found in the quasi-particle approximation. In the latter case we find a 
saturation point at the edge of the densities given, in the continuum method this is shifted 
to lower densities. The reason for this can be found in the next two columns, (K)/2 is 
the average kinetic energy, the first contribution to the energy in Eq. [2.20| . (V) /2 is the 
first moment of the spectral function (also known as the mean removal energy), the second 
term in Eq. |2.20| . In the continuum method the average kinetic energy is larger than in 
the quasi-particle method. This is obvious, in calculating (K) we integrate over n(p) with 
an additional weight p 2 over all momenta. Since strength is moved to higher momenta 
which have larger weights and (K) is larger than found in the quasi-particle approximation. 
Moreover, the amount of particles with a momentum larger than pf (n + (p)) slightly decreases 
with increasing density. Altogether this leads to an almost constant difference between 
the values of (K) calculated in the continuum and quasi-particle approximation. For the 
correlation contribution (V) the difference larger at the lower end of the calculated density 
range than at the higher end. Both effects combine to an effective repulsion which in the 
continuum approach shifts the saturation point to lower densities (see also Fig. |). We 
have to be careful in generalizing this finding: the relatively large shift in the saturation 
point might be a result of the system being only slightly bound which makes the saturation 
point very sensitive. Using a realistic interaction the binding will become stronger and the 
saturation point will be less sensitive to small changes in the self-energy. On the other hand, 
the stronger tensor force of such an interaction leads to a higher depletion and therefore 
enhances effects of the continuum contributions. 

In the fifth column we give values for the single particle energy at the Fermi-surface, 
equal to the chemical potential /i. In a thermo-dynamically consistent model this is related 
to the binding energy via the Hugenholtz-van Hove theorem (rjj : 

e sp (p f ) =E b + ^, (3.1) 
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where P is the pressure of the system, by definition P = p 2 dEb/dp. An obvious special 
case is that at saturation the chemical potential should be equal to the binding energy The 
difference between these two is a good measure of the thermo-dynamical consistency of the 
model. As we already noted above, the T-matrix approximation is a thermo-dynamically 
consistent approximation |1J. Since self-consistency is a crucial requirement for the thermo- 
dynamical consistency and the improved degree of self-consistency over the quasi-particle 
approximation one might expect improvements also with regard to the Hugenholtz-van Hove 
theorem. However, this is not the case, we find similar violations in both approaches: a 
violation at saturation of 5 MeV for both calculations. This is somewhat larger than we 
found for a relativistic calculation in the quasi-particle approximation [§],||, but in the same 
range. Although it is still far from ideal, it is already much better than the typical values 
found in the Brueckner approximation (i 20 MeV). 

In the two rightmost columns we give values for the density calculated using Eq. |2.21 



divided by the density calculated in the quasi-particle approximation (g^P/), divided in con- 
tributions with momenta below (n~(p)) and momenta above {n + (p)) the Fermi- momentum. 
The momentum integration was carried out up to 3pj, there n(p) is already much smaller 
than the accuracy with which we fullfill the spectral function sum-rule Eq. ^.8| . Particle 
conservation requires that the sum of the two contributions should be equal to unity. In our 
case the particle sum-rule is slightly too high, by a couple of promilles for both calculations 
with no significant difference between the two. Although this deviation is within our nu- 
merical accuracy there seems to be a trend that at higher Fermi-momenta, the calculated 
relative density is lower. A similar trend is also reported in the calculations of Refs [|3],[16 . 



Given the crucial role of self-consistency in either thermodynamical consistency and con- 
servation laws we find it surprising that we do not find significant differences between the 
violations in both the quasi-particle approximation and the continuum approach. Regarding 
self-consistency the continuum approach clearly goes beyond the quasi-particle approxima- 
tion. A fully self-consistent calculation would include a spreading of the quasi-particle 
peak. Since the widths involved there are not too large (most of the particles are around 



the Fermi-surface where the width is small), we think the effect of this would be smaller 
than the continuum contributions we included. Presuming this, it seems that (the lack of) 
self-consistency is not the reason for the violations we find. We already touched upon one 
possible reason, the validity of the gradient expansion. In Ref. |10] the conservation laws 
are proven in momentum space only within the gradient expansion, and one might speculate 
that the violations we find are caused by higher order derivatives of the T-matrix. We will 
investigate this point and possible improvements in the future. 

IV. CONCLUSIONS 

We presented a calculation for nuclear matter in which we went beyond the usual quasi- 
particle approximation and included off-shell effects in the self-consistent determination of 
the self-energies. We achieved this by a suitable approximation of the spectral function 
that preserves important properties like the sum-rule and the zero at an energy equal the 
chemical potential. For momenta below the Fermi-momentum (holes) our approximation 
amounted to putting all hole-strength into the quasi-particle peak and retaining the tail of 
the spectral function with particle strength. For particles a similar procedure was followed. 
It resulted in a model in which the additional continuum contributions are a perturbation 
on the quasi-particle contributions: in the limit of zero imaginary part of the self-energy the 
continuum contributions vanish and the model reduces to the conventional quasi-particle 
approximation. 

We presented results using a non-relativistic version of the Bonn potential in which 
we reduce the tensor force to avoid bound states. From the model calculation we found 
that the quasi-particle spectrum changes a few MeV upon incorporating the continuum 
contributions, the most for small momenta. The average strength typically present in the 
tail of the spectral function we now include is in the order of 6 %. In the equation of state we 
found that the continuum contributions provide an additional repulsion: the saturation point 
is shifted to a lower density. We also checked particle conservation and thermo-dynamical 
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consistency as expressed by the Hugenholtz-van Hove theorem. Since self-consistency is of 
prime importance in the proofs of these one would a-priori expect that the violations of 
these should diminish in the continuum approach because of the increase in self-consistency 
over the quasi-particle approximation. However, we found very similar violations for both 
approaches. We indicated that these might not be caused by lack of self-consistency but by 
higher order derivatives of the T-matrix, which are ignored in the gradient expansion. A 
better understanding and possible improvements will be subject of future studies. Another 
point of future attention must be the inclusion of bound states in the approach. This will 
allow to perform a calculation with a realistic interaction, these typically generate a depletion 
of ~ 15 % and one expects the effects of the continuum contributions to be even larger. 



APPENDIX: CALCULATION OF THE SELF-ENERGIES 

All integrations are transformed to cm. defining variables P and p. The incoming 
momentum is {u^i ab ,pi ab ), the momentum over which is integrated is (uJk,k), so that P = 
Piab + k,p= (pi ab — k)/2. The first contribution to the self-energy in Eq. ^4] (the Brueckner 
type term) becomes: 

T {+) g < (u labl p lab ) = 

[ P ™ x dP 2 I*™* dp 2 —^- (p\T^\u lab + uo Kqp )\p)Z h {k) + 

L dP d P (9^ / ^(p|T (+) (^ + ^)b)a(^,A;). (Al) 

JP Ln J Prm n {2ir) Z 2p lab J-OO 2lX 

The first integral integrates over the momenta k < pf, there the spectral function has a 
quasi-particle pole for Uk = ujk, qp - This is reflected in the limits of the integration: P^ in = 
mzx(p lab - p f ,0),P^ ax = p lab + p f and p 2 min<max = \{k 2 m ^ max + p 2 lab ) - ^ with k min = 
max(P — pi ab , 0), k max = min(p/, P + pi ab ). The second integral integrates over the momenta 
k > pf, there we have to integrate the energy variable over the tail of the spectral function. 
Here the integration limits are P^ in = maxfjo/- — pi ab , 0), P^ax = °°i the expressions for 
Pmin,max are tne same as before but now k min = max(P - p lab , 0),p f ), k max = oo. In the 
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practical calculation we limit this integration by k max — 2pf which implies P ma x = Piab+k ma x- 
In the calculation of Im\T < g^] we can identify various contributions. First we integrate 
separately over hole and particle external momenta. Each of these integrations have a quasi- 
particle part and a continuum contribution. Also T < has a quasi-particle and a continuum 
part. This all leads to eight different contributions to the self-energy. Note also that all 
integrations are restricted by the Pauli-blocking operator Qhh present in T < , specifically 
these are Q hh {P,p) = OVP > 2p f and Q hh {P,p) = OVp 2 > p) - P 2 /4. 

The most important contribution is the quasi-particle contribution from an integration 
over an external hole momentum with the quasi-particle part from T < , we denote this by a 
subscript qp,qp,h 

f Pmax f'Pmax I 

. dP 2 / dp 2 —— 

IP mi n J Pmin [2ltyApi ab 

(A2) 



. , P Pmax CPmax \ 

Im[T < g { ->(uji ab ,pi ab )} qPiqp , h = dP / dp (p\T<(u lab + u k<qp )\p)Z h (k). 



J min,max 



The integration limits are now P min = max(p /ab -p f , 0), P max = mm(p lab +p f , 2p f ), p Tl 
are the same as above and k min = max(P — p lab , 0), k max = min(p^, P + pi ab )- Recall that we 
also have the condition ui ab + LOk, qp = £^ + — + 2E(P, po) qp . Since the values which p can 
assume are limited by Qhh the region where the integrand is finite is further limited. We 
implement this numerically. 

We also have the equivalent contribution where we integrate over an external particle 
momentum 

, . C Pmax CPmax j 

Im[T < g { ->(ui ab ,pi ab )] qP:qp , p = dP / dp (p\T<(u lab + uj Kqp ) \p}Z p (k), 

J Pmin ^ Pmin 

(A3) 

for which we have the integration limits P m %n = max(pj — p lab , 0), P ma x — %Pf and k min = 
max(P — pi ab , 0),pf), k max = oo. Again in the practical calculation we limit this integration 
by k max = 2pf which implies P max = min(2pf,pi ab + k max ). The integration region is further 
limited by the limitations on p . From Eq. |2.19| it is clear that T < is zero for Q > 2/i. This 



results in the contribution above being zero for u>i ab = H since for particles u)k, qp > 
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A second class are the contributions where the external integration is over a tail of a 
spectral function with the quasi-particle contribution of T< 

Jm[T < 5( ( ~ ) (^ /afe ,p/, 



ab ) \qp,cont,h 

dP 2 I dp 2 —^- / ■ ~~ ^(p\T<(uj lah + u k )\p)a(uj,k), (A4) 



Pmax rPmax \ r2fl-Uli ab 

dP / dp TT^Ta / ~0 \P\ x qp 

Pm,„ ^Pmm (2lT) Z 4pi ab Jfl 2l{ 



1 mm 



and the integration limits are similar as for the qp,qp,h contribution, Eq. |A2| , P n 
max(pi ab - Pf,Q),P m ax = min(j9/ a b +pf,2pf), p m in,max are the same as above and k r , 
max(P — pi ab , 0), k max = min(pj, P + pi ab ). From the energy integration limits we see that 
this contribution is zero for Ui ab > fi. The particle momentum counterpart is: 

Im[T < g^\ui ab ,pi 



ab) \qp,cont,p 

dP 2 I dp 2 -—^- /' ^(p\T<(u lab + u k )\p)a(u,k). (A5) 



Pmax rpmax \ r Pmax dlO^ 

dP / dp 7T~vTa / 

With the same integration limits as the qp,qp,p contribution Eq. [Al| and ^ ma x = min(2/i — 

The third category are the contributions where we have the continuum contribution of 
T < with an external quasi-particle part 

, . C Pmax CPrnax \ 

Im[T < g { ->(uj lab ,pi ab )] cont:qP:h = dP / dp (p\T< nt (u lab + u km )\p)Z h (k). 

(A6) 

The integration limits are P min = max(pi ab -p f , 0), P max = pi ab +Pf, p m in,max are the same as 
above and k m i n = max(P — pi ab , 0), k max = min(p/, P + pi a b)- Its counterpart with a particle 
momentum quasi-particle contribution is 

, . rPmax rPmax \ 

Im[T < g^>(u lab ,pi ab )} cont ^ h = dP / dp (p\T< nt {u lab + u Kqp )\p)Z p {k). 

(A7) 

In the practical calculation we have the integration limits P m in = max(j>/ — pi ab , 0), P max = 
Piab + k max ) and k min = max(P - pi ab , 0),p/), k max = 2p f . Using the same argument as for 
the contribution Eq. |A3| it follows that this contribution is also zero for uji ab = fi. 
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In the fourth and final category we have the continuum contributions from both T< and 
the external integration: 

Jm[T < 5( (_) (tUlab, plab)]cont,cont,h = 

dP 2 / dp 2 —- / ^( p \T<(uj lab + u k )\p)a(u, k). (A8) 



Pmin Pmin 

' p, ZTT 

The integration limits are the same as in Eq. |A6| . From the integration limits it again follows 
that this contribution is zero for uii ab = \x. The particle-momentum counterpart is: 



Im[T < g { >{uJlab,Plab)]cont,cont,p = 

Pmax _ CPrnax _ j f f^rnax (],!<] 7, 

dP 2 / dp 2 —— / -±(p\T<(uj lab + u k )\p)a(u, k). (A9) 



'Pmin J Pmin {^) 2 ^Plab J ~OC ^Ti 9P 

The integration limits are the same as in Eq. [AT] and fi ma x = min(2/i — uJi ab ,n). 

In Eq. |2.9| we already presented the important condition that the imaginary part of the 
self-energy is zero at the chemical potential \i. In our calculation the imaginary part of the 
Brueckner type contribution Eq. [AT] is canceled by the various imaginary parts of T < g(~\ 
We can detail this further, we already identified the contributions to Im[T < g(~>] which are 
zero at uji ab = /x Using the unitarity relation of we can identify four contributions 
to the imaginary part of the Brueckner-type contribution. The first contribution of Eq. 



AJJ has a quasi-particle part integrated over a hole momentum, this is canceled by Eq. 



A2| , the continuum part integrated over a hole momentum is canceled by Eq. |A6| . The 
integration over a particle line of the quasi-particle contribution is canceled by Eq. |A5| and 
the continuum part is canceled by Eq. |A9|. 
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FIGURES 

FIG. 1. The quasi-particle spectrum calculated at normal nuclear matter density in the two 
approximations described in the text. The real and imaginary part of the self-energies at the 
quasi-particle peak in the continuum approach are represented by the solid and dotted lines re- 
spectively. The dashed and dash-dotted line stand for the real and imaginary part calculated in 
the quasi-particle approximation. 

FIG. 2. The real (lower panel) and imaginary (upper panel) of the self-consistent self-energy 
calculated in the continuum approach as a function of energy and momentum. For reasons of 
visibility the points of view are chosen differently and the imaginary part was multiplied by -1. 
Note the zero of the imaginary part at the chemical potential and the structure this induces in the 
real part. 

FIG. 3. The energy dependence of the self-energy for momentum equal to the Fermi momentum. 
The solid line is the real part of the self-energy in the continuum approach, the dashed line the 
real part of the self-energy in the quasi-particle approximation. The dotted and dash-dotted line 
are the imaginary parts of the self-energy calculated in the continuum and in the quasi-particle 
approximation respectively. On the scale of the figure the latter are almost identical. On the 
horizontal axis we use the variable u>', defined as the distance to the quasi-particle peak. 

FIG. 4. The binding energy and its components calculated in the two approximations discussed 
in the text. The solid line is the binding energy in the continuum approximation, the dashed line 
the binding energy in the quasi-particle approximation. The dotted and dash-dotted line stand 
for the kinetic contributions in the continuum and quasi-particle approximation respectively. The 
correlation contributions in the continuum and quasi-particle approximation are represented by 
the dot-dot-dash and the short dashed lines respectively. 
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TABLES 



Pf 


E b 


m 

2 


2 




n (p) 


n + (p) 


(GeV/c) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(-) 


(-) 


.25 


-4.2 


10.0 


-13.2 


-12.5 


0.935 


0.068 


.27 


-4.8 


11.7 


-16.5 


-12.6 


0.942 


0.061 


.29 


-5.2 


13.5 


-18.7 


-11.7 


0.946 


0.057 


.31 


-5.3 


15.4 


-20.6 


-10.0 


0.947 


0.054 


TABLE I. 


Various observables as explained 


in the text calculated in 


the quasi-particle approx- 


imation. 














Pf 


E b 


(K) 
2 


2 




n (p) 


n + {p) 


(GeV/c) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(-) 


(-) 


.25 


-3.4 


11.9 


-15.3 


-10.7 


0.934 


0.069 


.27 


-3.8 


13.6 


-17.4 


-10.5 


0.942 


0.060 


.29 


-4.1 


15.5 


-19.6 


- 9.5 


0.946 


0.056 


.31 


-3.9 


17.4 


-21.3 


- 7.3 


0.947 


0.054 



TABLE II. Various observables as explained in the text calculated in the continuum approxi- 
mation. 
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